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Abstract The light-by-light contribution from the lightest 
neutral pseudoscalar and scalar mesons to the anomalous 
magnetic moment of muon is calculated in the framework 
of the nonlocal SU (3) x SU (3) quark model. The model is 
based on chirally symmetric four-quark interaction of the 
Nambu-Jona-Lasinio type and Kobayashi-Maskawa-t'Hooft 
Ua(1) breaking six-quark interaction. Full kinematic depen- 
dence of vertices with off-shell mesons and photons in in- 
termediate states in the light-by-light scattering amplitude is 
taken into account. The small positive contributions from the 
scalar mesons stabilize the total result with respect to change 
of model parameters and reduces to aJ; bL ' PS+s = (6.25 ± 

0.83) -i(r 10 . 



1 Introduction 

The description of the muon anomalous magnetic moment 
(AMM) is one of the most challenging problem of the ele- 
mentary particle physics. Recent precise results on the muon 
AMM obtained in the experiment E821 at BNL [1] open 
possibility for very fine investigation of the contributions 
from the electromagnetic, weak and strong sectors of the 
standard model. At present, the theoretical predictions, based 
on e + <?~ annihilation and T decay inclusive cross sections, 
underestimate the experimental result by approximately 3 a 
(see, e.g. [2-4]). 

The most problematic part of the theoretical estimates is 
the contribution of the light-by-light (LbL) scattering thro- 
ugh the hadronic vacuum. The LbL scattering contribution 
(unlike to the hadronic vacuum polarization contribution) 
can not be extracted from phenomenological considerations 
or calculated from first principals. Different models were 



used for calculation of the LbL contribution to the muon 
AMM (see, e.g., [5-20]). 

In [16, 17] the contribution of the diagrams with pseu- 
doscalar meson exchanges to the muon AMM was estimated 
in the framework of nonlocal chiral quark model (N^QM) 
(for details see, e.g., [21-25] and Appendices of the present 
work). In the present paper the contribution of scalar meson 
exchanges is considered. 

Section 2 contains the description of the nonlocal chiral 
quark model, the meson dynamics in the pseudoscalar and 
scalar channels with mixing scheme and interaction with ex- 
ternal gauge field (see also Appendices A and B). The calcu- 
lation of the LbL contribution to the muon AMM from the 
pseudoscalar and scalar exchanges are detailed in Section 3 
(see also Appendices C, D and E). Our conclusions are given 
in Section 4. 



2 N^QM Lagrangian, T-matrix and r\ — r/' mixing 

The Lagrangian of the nonlocal SI/ (3) x SI/ (3) model is 



if = q(x)(id - m c )q(x) + - [J a s (x)J a s (x) + J a P {x)J a P {x) 



H 

-£ T abc 



J§(x)J b s (x)J$(x) - 3r s {x)J b P {x)J P {x) 



(1) 



where q (x) is the quark field, m c is the diagonal matrix of 
the quark current masses 1 , G and H are the four- and six- 
quark coupling constants. Last line in the Lagrangian repre- 
sents the Kobayashi-Maskawa-t'Hooft determinant vertex 
[26, 27] with the structural constant 
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(2) 



1 We consider the isospin limit m c u = m c j ^ m c s . 
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where X a are the Gell-Mann flavor matrices for a = 1,..,8 
and Ao = ^2/31. The nonlocal quark currents are 

J a ch {x) = Jd 4 x l d 4 X2f{x l )f(x 2 )q(x-x l )r^q(x + x 2 ), (3) 

where ch = S,P and IT? = X" for the scalar channel, = 
iy 5 X" for the pseudoscalar channel, and f(x) is a form factor 
reflecting the nonlocal properties of the QCD vacuum as it 
occurs in the instanton liquid model. 

The model can be bosonized using the stationary phase 
approximation which leads to the system of gap equations 
for the dynamical quark masses m d (i = u,d 7 s) 2 



m d ,u + G5„ + — S U S S = 0, 



m d ,s + GS S + —Si = 0, 



Si = -8AT C 



d%K f 2 (K 2 )m,(K 2 ) 
(2nr D,(K 2 ) ■ 



(4) 



where m^K 2 ) = m c j + m dti f 2 (K 2 ),Di(K 2 ) = K 2 +mf(K 2 ), 
f(K 2 ) is the nonlocal form factor in the momentum repre- 
sentation. 

The vertex functions and the meson masses can be found 
from the Bethe-Salpeter equation. For the separable interac- 
tion, given by Eqs. (1), (3), the quark-antiquark scattering 
T-matrix in the pseudoscalar (scalar) channel becomes 



T cA = f ch (P 2 )S 4 (Pi +P2-P3-P4) l\f( p ?)> 
1 



4 

n 

i=l 



-G^+n ch (P 2 



2\ I ^ ch ' 



(5) 



where P, are the momenta of external quark lines, G c /, and 
n c /,(P 2 ) are the corresponding matrices of the four-quark 
coupling constants and the polarization operators of mesons 
(P = P\ +P2 = P3 +P4). The meson masses Mm are deter- 
mined from the zeros of determinant, det^Gjl 1 — n c /,(— M^)) : 
0. The actual expressions for the matrices G c /, and Yl^ are 
given in Appendix A. 

The T-matrix for the system of mesons 3 in each neutral 
channel can be expressed as 



t cIi (p 2 ) = Y, Vm V 



M 



-(P 2 - 



(6) 



2 Through the paper the capital letters will be used for Euclidean mo- 
menta, while small letters for Minkowski momenta. 
3 Such description of the light scalar mesons as (^/-states is probably 
simplified. It seems that it is necessary to include other structures, e.g., 
four-quark states (see, e.g., [28]). However, the present model is for- 
mulated in the leading order of the 1 / N c expansion and our calcula- 
tions are consistent within given approximation. Moreover, the scalar 
mesons participate in the processes under consideration only as inter- 
mediate states, being far from mass-shell. 



where Mm are the meson masses, Vm{P 2 ) are the vertex 

functions (Vm(p 2 ) = "/ > V M (P 2 )'f l J . The sum in (6) is over 

full set of light mesons : (M = n° , 77 , 77 ' ) in the pseudoscalar 
channel and (M = ao(980), <7,/o(980)) in the scalar one. In 
general case of three unequal quark masses it is necessary to 
solve the n° — 77 — 77/ and qq — a — fo systems. However, in 
the isospin limit considered here they reduce to the n° and 
77 — 77' systems and to the ciq and a — fo systems. Then, it is 
convenient to diagonalize the scattering matrix by orthogo- 
nal transformations 



77 \ / cos Bp - sin dp \ ( 778 
777 \sin0p cos0/> ) \T\o 

a \ _ /cos d s - sin G s \ ( a 8 
/ (980y - ^sin^ cos0 5 J \a 

As a result the mesonic vertex functions are 



V ao (P 2 )=ig ao (P 2 )X 3 , 

V a (P 2 ) = ig a {P 2 ) (A 8 cos0 5 (P 2 ) - A()Sin0 5 (P 2 )) , 
V fo (P 2 ) = ig h {P 2 ) (A 8 sin e s (p 2 ) + A cos 9 S (P 2 )) , 



(7) 



(8) 



where gM(P 2 ) and 9 (P 2 ) are the meson renormalization con- 
stants and mixing angles (A. 5) depending on the meson vir- 
tuality. The renormalization constants are defined through 
the unrenormalized meson propagators Dm{P 2 ) as 



g 2 M (P 2 ) = -(P 2 +K)D M (P 2 ) 



(9) 



The meson mixing angles depend strongly on the meson 
virtuality. Therefore 9 a = 9 S (-M 2 a ) and 9 fo = O s (-M 2 fo ) 
are different for the on-shell a and fo mesons. The same 
situation takes place in the pseudoscalar sector, where 9^ = 

d P (-M 2 ), e n , = d P (-M 2 ,). 

External fields are introduced by derealization of the 
quark fields q(x) by using the Schwinger phase factor E(x,y) 



q(y) -> Q{x,y) = E(x,y)q(y), 

where 



(10) 



E(x,y) = ^exp i i / dz^(z) + ^{z)j 5 ]T a ), (11) 



and and are the external vector and axial-vector 
gauge fields, T" = X a /2. The ^exp is handled with help 
of prescription for the derivative of contour integral 

y y 
^ Jdz v F v (z)=F,(y), 8^ (x-y) J dz v F v (z)=0, 

X X 

as described in [29]. As a result the kinetic part leads to 
usual local electroweak vertices. However, the terms with 
nonlocal quark currents J" h (x) generate additional vertices, 
see Appendix B. 
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(a) (b) (c) 

Fig. 1 LbL contribution from intermediate meson exchanges. 

For numerical estimates we use the Gaussian nonlocal 
form factor for Euclidean momenta f(K 2 ) = exp(— K 2 /2A 2 ) 
and the model parameters obtained in [30]. The model pa- 
rameters (the current quark masses m c j, the coupling con- 
stants G and H, and the nonlocality scale A) are fixed in 
[30] by requirement that the model reproduces correctly the 
measured values [31] of the pion and kaon masses, the pion 
decay constant f n , and the 77' mass (parameter sets G/, Gjv) 
or the 77' — » 77 decay constant g^yy (sets G//, G///). The sets 
G/, Giv vary by different input for the nonstrange current 
quark mass, while G//, G/// are two solutions of the same 
fitting procedure. 

3 LbL contribution from resonance exchanges 

The basic element for calculations of the hadronic LbL con- 
tribution to the muon AMM is the fourth-rank light quark 
hadronic vacuum polarization tensor 

IWpfei.ft.ft) = / A, j d 4 x 2 j d 4 x 3 e' iinx>+!nX2+c ' 3X3] x 

(12) 



x(0\TUn(xi)jv(x 2 )h(x3)j P (0))\0) 



where (x) are light quark electromagnetic currents and |0) 
is the QCD vacuum state. 

The muon AMM can be extracted by using the projec- 
tion [32] 



where 



t— r / / \ . 6 [ ^ 4 gi f d4 q2 1 



x 7 



{p'-qi) 2 -ml {p-qi-q 2 ) 2 -m 2 l 



X -^-p^tivXaiqi, q2,k-qi-q 2 ), (13) 
'in is the muon mass, = (p' — p)^ and it is necessary to 



consider the limit kn — > 0. 



In the case of the resonance exchanges of the light hadrons 
in the intermediate pseudoscalar and scalar channel the LbL 
contribution to the muon AMM is shown in Fig. 1. The ver- 
tices containing the virtual meson with momentum p and 



two photons with momenta q\, 2 and the polarization vectors 
£12 (see Appendix C and Fig. 6) can be written as [15] 



at 

with 



(92^2) 



(j>) 



e 2 e?E 2 *A p lV (p,q h q 2 ) (14) 



A£ V (p,q h q 2 ) = -ie^ vpa q P q 2 F P (p 2 ;q 2 ,q 2 ) 
and 



(15) 



^?i,ei)Tf< 



(12, £2) 



V: 



e 2 £?£ 2 v Ag v {p,q h q 2 ) (16) 



with 



&s v (.PiixAi) = As (p 2 ;q 2 i,q 2 v )PA V (9i,q2)+ 

+B s (p 2 ;qlq 2 2 )Pg v (q u q 2 ), 



(17) 



where 



P A V (9hq2) = (g^ v (qiq2)-qiq 2 ), 

fs v (<7i,<72) = (q 2 q 2 - (9192)91) {q\q\ - {qm)ql) , 

and p = qi + q 2 - Note, that the scalar form factor B5 is sin- 
gular in the limit when one photon is real and the virtuality 
of the second photon equals to the virtuality of the scalar 
meson p 2 — > q\, q\ — ^ 0. For convenience we also define an 
additional function 

B' s (p 2 -qlql) = B s (p 2 -q 2 u q 2 ) {{q x q 2 ) 2 - qjqj) , (18) 

which is regular in this limit. In general case, these scalar 
functions are combinations of the nonstrange and strange 
components. Details for the mixing of mesons interacting 
with photons are given in Appendix C. In Appendix D the 
local limit of the amplitude 7*7* — > S* is presented. 

The polarization tensor Tl^^P for the exchange of me- 
son with mass Mm is 



Yl^P( qi ,q 2l q 3 ) 



(19) 





^q 2 ,q u q 2 )A l P(qi + 


<72,<73,<74) 




(?i+ft) 2 -M M 2 




AW(q 2 - 


>rq3,q\Ai) AVX ("724 


q3,qi,q3) 




(q 2 +q 3 ) 2 -M M 2 




A^(qi- 


f43,<7i,<73M vp (<7i -+ 


q^AiA*) 


{qi+q3) 2 -M M 2 



+ i 



where qi are momenta of outgoing photons, ^4 = — (q\ + 
q 2 +qi), and one should take Ap for pseudoscalar and A £ 
for scalar mesons, respectively. Details for the pseudoscalar 
exchange can be found in [10] (Eqs. (3.1) and (3.3)). The 



4 



low-energy expansion of the derivative of the polarization 
tensor Y\P- V ^P for the scalar meson is given by 



dkP 



Yl' lvX(J {q u q 2 ,k-q l -q 2 ) = 



Ag V (q l +q 2 ,qi,q2) d ^ A 
1 ( qi +q 2 ) 2 -M M 2 dkP S 

As X (-n,i2,-qi -qi) d 
dkP 



h +q2,-qi-q2>-k) 



q\ - M M 2 



Ag"(-q u q u -k) (20) 



+ i 



As X (-q2,qi,-qi-q2) d 



q z 2 -M 



M 



dkP 



A™(-q 2 ,q 2 -k) + 0(k). 



and the low-energy expansion for the derivative of A^ v is 



d 

dkP^ s 



4f (-q, q, k) = A s (q 2 , q 2 , 0) ( g» v q» - q v g^ 
q^q 



B' s (q 2 ,q 2 ,0)q" 



r 



0(k). 



(21) 



As a result the numerator of the two-loop integrand for the 
aJ; bL contains the combination of two form-factors and is a 
polynomial in momenta. 

At next step, the expression for LbL can be averaged [2] 
over directions of the muon momentum p 

, , 1 



2tt 2 



dQ (p) . 



(22) 



After averaging the expression for the LbL contribution to 
the muon AMM from the light scalar meson exchange can 
be written in the form of integral over Euclidean momenta 

4 bL ' S - ]dQ\ JdQl J dtVT^it: 
o o -1 



Qi 



(23) 



<< = E I Xs(QhQiO)Y s (QiQ 2 ,Q 2 )T4\ 

X=A,B'Y=A,B 

= I I Xs(QhQi0)Y s (Q 2 ;Q 2 uQ 2 )Tsr, 

X=A.B'Y=A,B 

where X$, Y$ are the functions A$, B$, B' s deined in Eqs. (17), 
(18), 23 = — (Qi + Q 2 ) and capital letters are introduced for 
Euclidean momenta, i.e. Q 2 = —qf. One should note that 
Tsf' A = iTsf A , Tsf' B = jTs^ B . The functions Tsf Y are 
given in Appendix E. 




150 



200 



250 300 

m,i, u [MeV] 
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400 



Fig. 2 LbL contribution to the muon AMM from the neutral pion 
and o" exchanges as a function of the dynamical quark mass. Bunch of 
three lower lines correspond to the o" contribution, the 71° contribution 
is in the middle, and the upper lines are the combined contribution. The 
band along the thick line between dashed and dotted lines corresponds 
to the error interval for the pion two-photon width. Vertical thin dashed 
lines denote the interval of dynamical quark masses used for the esti- 
mate of the error band for «J; bL . 



where the corresponding flavor matrices are A u = (V2Xq + 
^)/v / 3=diag(l,l,0) wdb= 1,2,3. In this case the model 
has three parameters: the current quark mass m c u , the dy- 
namical quark mass rrid u and the nonlocality parameter A . 
In order to understand the stability of the model predictions 
with respect to changes of the model parameters one may 
vary one parameter in rather wide physically acceptable in- 
terval, while fix other parameters by using as input the pion 
mass and the two-photon decay constant of the neutral pion. 
Thus, we take the values of the dynamical quark mass in the 
typical interval of model values 200-350 MeV and other pa- 
rameters are fitted by the above physical observables within 
the error range given in [31]. 

The results are shown in Fig. 2. We see that the contri- 
bution of the a meson is small, positive and has very small 
minimum value around the value of 300 MeV for the dy- 
namical mass. It is interesting to note that the total result 
for the pion and a meson contributions is rather stable to 
variation of the dynamical mass in the tested interval. Our 
estimates for the n° and the sum of n° and a contributions 
(here and below in 10~ 10 ) are 



LbL,7T° 



= 5.01 ±0.37, 



LbL,!r°+(T 



= 5.40±0.33. 



(24) 



4 The results of model calculation 

It is instructive to study the n and a mesons contribution 
to the muon AMM for the SU(2) x 51/(2) version of the 
nonlocal model. The Lagrangian of the model is given by 



: q(x)(id - m c )q{x) + | [j u s (x)J" S (x) + 4(x)J b P (x) 



In the 51/(3) model, for the central values of 77 and T]' 
contributions we use the averages over different parameter- 
izations [17]. The error bar for r\' is taken as a maximal 
deviation from the central value. The deviation of the 77 con- 
tribution from the central value seems accidentally small, so 
we use the factor 60% of the error value for tj' as an estimate 
of the error bar for the tj contribution. Our estimate for the 
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set 


7t° 


n 


1' 


rj+f}' 


PS 


«o(980) 


O" 


/o(980) 


S 


x° + o 


PS+S | 


G, 


5.05 


0.55 


0.27 


0.82 


5.87 


0.0064 


0.100 


0.0035 


0.110 


5.15 


5.98 


Git 


5.05 


0.59 


0.48 


1.08 


6.13 


0.0079 


0.100 


0.0038 


0.110 


5.15 


6.24 


Gin 


5.05 


0.53 


0.18 


0.71 


5.76 


0.0058 


0.100 


0.0034 


0.109 


5.15 


5.87 


G/v 


5.10 


0.49 


0.25 


0.74 


5.84 


0.0060 


0.115 


0.0038 


0.126 


5.25 


5.97 



Table 1 The contribution of scalar and pseudoscalar mesons to the muon AMM ajf for different sets of model parameters [30]. All numbers are 
given in 10 . 



77 and rj' contributions is 

a^ L ' n =0.54 ±0.32, a^' 11 ' = 0.30±0.18. (25) 

Finally, we estimate the combined contribution from the 
a (980) and / (980) mesons as 



aJ^ L - ao+fe w0.01, 



and add it to the total result. 



(26) 



5 Comparison with other models 

It should be mentioned that there are estimates of the scalar 
meson exchange contributions to the muon AMM in differ- 
ent versions of the local NJL model. In [ 14] the combined 
scalar contribution was estimated as 



a 



„LbL,S 



(0.68±0.2)-10 



-10 



(27) 



whereas [15] gives the estimates for a and «o contributions 

as 



a 



a 



LbL.cr 
M 

LbL.ao 



= (1.167±0.238)-10 
(0.062±0.024)-10 



-10 



-10 



(28) 



One can see that our results (Table 1) are smaller in absolute 
values than other estimates. 

Note, that in estimations (27) and (28) there is an am- 
biguity in the sign for the scalar meson exchange contribu- 
tions. In [33] the analytical expressions for the pion and a- 
meson contributions 4 was obtained with the meson transi- 
tion form factors taken from the simple vector meson dom- 
inance (VMD) model parameterized by the p -meson mass 
M p . These expressions are given as an expansion in small 
parameters chosen in accordance with the mass scale hierar- 
chy (My — rnjj ) <C m 2 ^ <C M 2 , . We reproduce numerically the 
coefficients of these expansions for the pseudoscalar meson 
contribution (Eqs.(8) and (10) of [33]) 5 by using our code 



4 For the pion exchange contribution, the coefficient of the leading, 
log 2 (Mn/mu), term in (m^/Mp) 2 expansion was found in [34]. 
5 Namely, varying the mass parameter values for the p -meson, muon 
and difference between the pion and muon masses squared, one can 
extract different terms of the expansions given in [33]. 



for numerical calculations of and substituting the VMD 
transition form factors instead of NjQM ones 



( 2 2 2\ 



- kVMD I 2 2\ 



Bs 



( 2 2 2\ 
[P -1ul2) 



0. 



However, for the c-meson contribution by using the same 
VMD model as in [33], we reproduce numerically the coef- 
ficients of the expansion in Eq. (1 1) of [33] up to the overall 
sign. Thus we conclude (as in [15]) that the scalar meson 
contribution to the muon AMM for the VMD model has a 
positive sign in variance with the result of [33]. 

For additional check of computer code correctness, as 
suggested in [33], one can calculate the contribution to the 
muon AMM from the vacuum polarization processes, a^ 1VM ^ , 
where the virtual photon splits to the meson M and the real 
photon y. As argued in [33] these contributions have to be 
positive, since they are related by dispersion relations to the 
cross sections o(e + e~ — > My). Our numerical results for 
the VMD model and a set of parameters used in [33] (for 
M a = M n ) are 



a^£ = +5.64.10- 10 , 



>,VMD 

LbL,CT 
2 11, VMD 
_0, 



n = +4.76-10- 10 , 



VvMD r = +0.368- 10- 10 , 
VvMD r = +0-265. 10-">. 



(29) 
(30) 

(31) 
(32) 



Our numerical results for the pion exchange contributions, 
(29) and (31), agree with the results given in Eqs. (2) and 
(10) of [33]. 

In order to study the transition from the nonlocal model 
to the local one, we consider 6 the nonlocal SU (2) x SI/ (2) 
model with Pauli-Villars regularization parameterized by 

1 . parameter of nonlocality A , 

2. parameter of quark loop regularization A q . 



The local model corresponds to the limit 



A 



A,=A, 



fit 



(33) 



while the nonlocal model without regularization can be ob- 
tained by setting 



A=A m , A„ -> ■ 



(34) 



Similar consideration was used in [25] for the investigation of l/N c 
corrections. 
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9.0 




ffl — 

3 5.0 
% 4.0 
g 3.0 

2.0 \ 

1.0 

o.o I — 1 — 1 — 1 — " r 7 ' 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

x 

Fig. 3 LbL contribution to the muon AMM from the neutral pion and 
(7 exchanges in nonlocal model with additional regularization. Zero x 
corresponds to a local NJL model while x equal 1 is the nonlocal model 
without regularization. 

For definiteness, let us compare 7 the local model 8 [35, 36] 
with the nonlocal one 9 [37]. The values of the quark conden- 
sate in the local and the nonlocal models coincide numeri- 
cally within less than 0.5 % deviation. 

The result is presented in Fig. 3, where we introduce the 
parameter of nonlocality x. Zero of x corresponds to the local 
NJL model, while the nonlocal model is reproduced for x 
equal to one. For each x point in between zero and one the 
regularization parameter behaves as A q = A?° L /(1 —x), the 
dynamical quark mass scales linearly and the current quark 
mass and A are refitted in order to reproduce the mass of the 
neutral pion and linearly scaled quark condensate. 

One can see that the pion contribution to the muon AMM 
increases from 5.3 • 1CT 10 in the NjOM model to 8.5 ■ 1CT 10 
in the local NJL model. More dramatic situation takes place 
for the <7-meson contribution. The contribution in the local 
limit is ten times larger than in the nonlocal model (2.2 • 
10" 10 instead of 0.22 ■ 10" 10 ). The values obtained in the 
local limit are of the same order as numbers quoted in [ 15]. 



6 Conclusions 

We found that within the N^QM the pseudoscalar meson 
contributions to muon AMM are systematically lower then 
the results obtained in the other works (see discussion in 
[17]). The full kinematic dependence of the vertices on the 
pion virtuality diminishes the result by about 20-30% as 
compared to the case where this dependence is neglected. 
For 77 and 77' mesons the results are reduced by factor about 
3 in comparison with the results obtained in other models 



We rescale current quark masses in order to reproduce mass of neutral 
pion instead of charged one. 

8 m™, L = 5.69 MeV, = 253.9 MeV, A ? NJL = 800 MeV. 
Vi yV = 5.45 MeV, m** QM = 255.8 MeV, A N ^ M = 902.4 MeV. 



where the kinematic dependence was neglected (see details 
in [ 17]). The total contribution of pseudoscalar exchanges 

a^ L ' VS = (5.85 ± 0.87) • 10 -10 (35) 

is approximately by factor 1 .5 less than the most of previous 
estimates. 

The scalar mesons contribution is positive and partially 
cancels model dependence of the pseudoscalar contribution. 
The combined value for the scalar-pseudoscalar contribu- 
tion is estimated as 

a LbL,PS+s = (6.25 ± 0.83) -10 -10 . (36) 

It is well known that the hadronic LbL contribution to 
the muon AMM calculated in the effective approaches have 
model dependent features. However, there are few constraints 
following from QCD and need to be satisfied in any accept- 
able model calculations. One result was obtained by Mel- 
nikov and Vainshtein [11] (see also for discussions [38]) for 
the high photon momentum behavior of the total LbL am- 
plitude. The consistency of the N^QM with the Melnikov- 
Vainshtein constraints was carefully analyzed in our previ- 
ous works [16, 17]. It is based on the fact that in the high 
photon momentum limit the dynamical nonperturbative dress- 
ing effects containing in the quark box diagram becomes 
vanishing and this diagram produces short range behavior 
characteristic for the perturbative QCD regime. At the same 
time the LbL contributions containing the meson exchanges 
are responsible for the long distance dynamics and are sup- 
pressed in the high photon momentum limit. 

Another important constraint concerns the low photon 
momentum limit of the LbL contribution with intermedi- 
ate pion exchange. In [34] (see also for discussions [38]) 
the coefficient of the leading logarithm of ultraviolet regula- 
tor, log 2 A, arising in this contribution was computed in the 
chiral limit in the leading in 1 /N c approximation. In [34] 
there was also discussed the correspondence between this 
result and the model calculations given in [10]. They found 
that for the case of the vector meson dominance form fac- 
tor and in the limit My —¥ 00 the logarithmic coefficient nu- 
merically agrees with chiral perturbative theory result. We 
have checked this statement by performing the chiral expan- 
sion of the polarization operators as well as triangle diagram 
functions. Then, one can indeed reproduce the correct coef- 
ficients in front of log 2 A terms in the limit as the ultraviolet 
regulator A goes to infinity. 

Finally, the important point for the model calculations is 
the total contribution from all leading diagrams. This is be- 
cause the different models may redistribute partial contribu- 
tions differently. In the present work we show that the con- 
tribution of the scalar mesons being relatively small leads 
to stabilization of the total pseudoscalar and scalar contribu- 
tion with respect to variation of the model parameters. The 
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next step is the calculation of the dynamical quark and pion 
loops contribution which is now in progress. 

We thank M. Buballa, Yu. M. Bystritskiy, C. Fischer, 
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channel. The unrenormalized mesonic propagators for the 
scalar mesons are 



D; 1 (/ ,2 ) = -G ao 1 +n ao (p 2 ), 



(A + C)±J(A-C) 2 + 4B 2 



A = -G Hii /det(G ch )+n ()0 {P 2 ), 
fi = +G 8/det(G c . /; )+n 8( J P 2 ), 
C=-G 00 /det(G ch )+Ihs(P 2 ), 
det(G c ft) = GooGgg — ^08- 

The mixing angle depends on the meson virtuality 



(A.4) 



Appendix A: Four-quark coupling constants, 
polarization operators and mixing angles 

The elements of G c /, -matrices take the form 



Goo = G ± % (2S U + S s ), G 88 = G T ? ( 45 » - s *) . 
3 6 

G m = G m = ^^-H(S u -S,), G 33 = Gt^S s , (A.l) 
6 I 



where the upper sign corresponds to the scalar channel, while 
the lower sign corresponds to the pseudoscalar channel. For 
the pion, G K equals to G33 of the pseudoscalar interaction, 
and, for the ao-meson, G aQ equals to G33 of the scalar inter- 
action. 

The elements of n c /,(P 2 ) -matrix for the scalar and pseu- 
doscalar mesons are diagonal in the quark-flavor basis, and 
in the singlet-triplet-octet basis they are given by 



n 0Q (p 2 ) = ^(2n uu (p 2 )+n ss (p 2 )), 

n 88 (p 2 ) = i(n !( „(p 2 ) +2n, s (p 2 )), (A.2) 

n 0H (p 2 ) = n 80 (p 2 ) = ^(ii„„(p 2 ) - n ss (p 2 )), 
n 3i (p 2 ) = n uu (p 2 ), 

where the difference between the scalar and pseudoscalar 
channels is in the polarization operators 



17,,-OP 2 ) =8/Y, 



4K f 2 {K\)f 2 {K 2 _) 



(2nyD,(K 2 )Dj(K 2 ) 
x [{K + -K-)^ mi (Kl)mj(K 2 )] , (A3) 

where K± = K±P/2. Similarly to Eq. (A.l), the upper sign 
corresponds to the scalar channel and the lower sign corre- 
sponds to the pseudoscalar channel, TI ao equals to 7733 for 
the scalar channel and Tl n equals to 1133 the pseudoscalar 



s (P 2 ) = iarctan 



IB 



A-C 



-&(A-C). 



(A.5) 



Expressions for the unrenormalized propagators for the 
pseudoscalar mesons are similar to the scalar meson propa- 
gators, Eqs. (A.4), (A.5), with replacements oq K> f ?7 . 
f -> ??' and 9 S -> dp. 



Appendix B: Feynman rules for nonlocal vertices 




k » 

Pi Vi 

(b) 



Fig. 4 Vertices r£ 2iPl , Eq. (B.6), V^ puq , Eq. (B.8), with one photon. 



The total vertex of photon interaction with quark-antiquark 
pair (Fig. 4a) contains local and nonlocal parts 



rL Pl = T a (Yn-(pi+p2)nT a W(p 1 , P2 )), 



(B.6) 



where T" = Q and m^\p\,p 2 ) is the first order finite-difference 
of the dynamical quark mass 



m W(p 1 ,p 2 ) = mi ^ ) - m(P2) 



2 2 

pi -pi 



(B.7) 



The contact interaction vertex of meson, photon and quark- 
antiquark pair (Fig. 4b) is purely nonlocal and takes the form 

rf 2 ; ^ =~gu{k) (^(pi,Pi + q)f( P 2)(2p 1 + q)> l T a r»+ 
+f^(p 2 ,p 2 -q)f( Pl )(2 P2 - q )^r') . (B.8) 
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92(f) 



Pi P2 

(a) 




Fig. 5 Vertices T P2 ^, uq ^ qi , Eq. (B.ll), and rj££V<», Et l- ( B - 12 )- 
with two photons. 



In order to express the vertices with two external photons 
we introduce the following functions 

G a ll (k,q)=iT a (2k + q)^(k,k + q), 
Gf v (k )q ,q',k') = -f(k') jrr 6 [g^ l \k,k + q + q ')+ 
+ [2(k + c/)+ q ]^(2k + q ') v fW (k,k + q ',k + q + c/)] + 
+ [(q,a, l l)<—>(q',b,v)]}, (B.9) 

where f( 2 ) {k\^k 2 ,k^) is the second order finite-difference 

,W MjW . gfe^M , (B.10, 

3 2 

With this notation, the vertex of two-photon interaction with 
quark-antiquark pair (Fig. 5a) is 



-m d {Gf v {pi, q i, q2 ,P2) 



+ Gf v (p2- q \- q 2, q \, q i,P\)+ (B.ll) 

+G£ Ot , q 1 ) G* (p 2 - <12, 92) + G£ (p 2 - q 1 , q 1 ) G\ (pi , q 2 ) 

and the interaction vertex for two photons, meson and quark- 
antiquark pair (Fig. 5b) becomes 

^Pimrn =-SM{k){Gf v ( Puquq2 , P2 )r^ 



+ r M G 1i V {p2-q\ -<i2,qi,q2,pi) 

+G a tl (pi, q i)r£G b v (p2- q 2,q2) 

+G^ip 2 - q i, q i)r^G b v (p 1 , q2 ) 



(B.12) 



meson it takes the form 

A(f l f 2 ^M*)=e 2 e^A^(p, quq2 ) 



A^ v (p,quq 2 ) = -iN c 



d 4 k 



(C.13) 



x Tr (21^ S(h )r£ k3 S(k 3 )Tl >k2 S(k 2 ) 

where the symbols are the photon momenta 171.2-. the photon 
polarization vectors E\ i2 , the meson momentum p = cji +cj 2 , 
and the quark momenta £1,2.3 (ki=k + q \,k 2 =k — q2 ,ks = 
k). The first term in parentheses corresponds to the quark 
triangle diagrams 10 (Fig. 6b and crossed term Fig. 6c) and 
next terms corresponds to the diagrams in Figs. 6d-g with 
effective nonlocal vertices defined in (B.8), (B. 1 1), (B.12). 

For different scalar meson states one has the following 
combinations of nonstrange and strange components 

^ V (P,qi,q2) = 8a ip 2 )Sr(p 2 ;qi,ql), 



K {p,qi,qi) 



gq(p 2 ) . 
3V3 



(58r(pW l ,4)-28r(p 2 ; q W2))cos9 s ( P 2 )- 

^(58r(p 2 iqi,4) + 8 s ^(p 2 ;qlq 2 2 ))^e S (p 2 ) 
8f (P 2 )„ 



(C14) 



(55,r tf;&ti)-28r {p 2 ;qUl)) ™e s (p 2 ) 

V2(5sr ( P 2 ; q w 2 )+8r (pW uq i))co S e s (p- 



One can easily see from Eqs. (17), (18) that the mixing 
for the form-factors As, B5, B^ from the components A u , 
B u , B' u and A s , B s , B' s is similar. One should project A, and 
Bi (i = u,s) from loops of nonstrange and strange quarks 



Ai(p 2 ;q 2 ,q 2 ) 



2{q\-qi) 



(31-92) ?T?2 



(31 -<72) 2 - q\q\ 



B i (p 2 ;qj,qi)=- 



$r (p 2 ,q 2 ,q 2 ) 



X 



2(91-92) ((91 -92) 2 -9i92) 
(91-92) 9[92 



(91 -qi) -q\q\ 



where 

\iv (Jl „2 „2\ 



Appendix C: Amplitude with meson and two photons 8f (p , q 2 ,(} 2 ) = — 2i 



<Pk 



J bc +J de +J f 



(C.15) 



- J g V 



The photon-meson transition amplitude is a sum of diagrams 
shown in Fig. 6, where all particles are virtual. For the scalar 



LO 



In the case of pseudoscalar mesons, the diagrams in Fig. 6d-g give a 
zero contribution due to chirality considerations. 
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(a) (b) (c) 

Fig. 6 The diagrams for the photon-meson transition. 




(e) 




(g) 



and different terms in brackets, J^ v , correspond to the dia- 
grams shown in Figs. 6, with lower indices being the symbol 
of the figure. 

Below for simplicity, a momentum is denoted as a lower 
index and a quark flavor index i is omitted 

/„ = f(kl), m„ = mt{i*), D n = fj& = f W (k,„k m ), 

= m^(k n ,k m ), Q = (k n ,k m ,k,) . 



'be 



Then, one has 

hh 

D X D 2 D 3 



Vr + (k 2 + k 3 ) v (k l + k 3 fV A - 



mg (** + ~ m l3 (*i + k lT V3 



(C.16) 



where 



V? v = mi [klk^+k^kl] + m 3 [kltf-k^X] + 

+ m 2 [k 3 k^ +£3 k\~\ + g^ v [m\m 3 m 2 — m\(k 2 k 3 )- 
+ m i (k 2 ki) -m 2 (k 3 ki)] , 
V 2 = fcf [m 2 m 3 + (k 2 k 3 )} +k% [ni\m 3 — (&1&3)] + 
+ k% [m\m 2 + (kik 2 )] , 

= k\ [m 2 m 3 — (£2^3)] +k 2 [m\m 3 + (k\k 3 )] + 
+ ki[m 1 m 2 + (k l k 2 )}, 
V4 = m\m 2 m 3 +m\ (k 2 k 3 ) + tn 3 (k\k 2 ) + tn 2 (k\k 3 ), 



and 



x (rn 2 k v 3 + m 3 k v 2 -m$(k 2 + k 3 ) v ((k 3 k 2 ) + m 3 m 2 ^j + 



(C17) 



x ^m 3 ^ + m\k% —va^(k\ +k 3 )^((k 3 ki) + m 3 m\)^ , 



hh 

D X D 2 



{{k\k 2 ) + tni?n 2 )tn c i 



(/l+/2)^ V f ( ,2 ) + 



+ (*, +k 3 r(k 2 +k 3 y ((/, +/ 2 )(f 3 ) 1 +f^)-f®f$) 



r = -^rt 5 + (*i +hY (k2+hy^ l ] 



D 2 



l^ [g ^ + (k 2 +k 3 y(k l+ k 3 y^ 2i 



'12 

+ i^ifWfW {2 k l+ k 3 y(2k 2 +k 3 y 



D 3 



Analytical expressions for the form factors A,- {p 2 ; q\ , q 2 ,) 
and B\ (p 2 ',qi,q 2 ) in the case of special kinematics 11 , when 
one photon is real, q\ = 0, and the virtuality of second pho- 
ton is equal to the virtuality of meson p\ = q\, can be ob- 
tained by expanding the quark-loop expressions, Eqs. (C. 15), 
(C.16), (C.17), in q\. The resulting expressions contain deriva- 
tives of the nonlocal function f(k 2 ) up to third order. These 
expressions are rather cumbersome and not presented here. 
Alternatively, one can calculate the form factors for small 
but nonzero q\ and then take the limit numerically. 



Appendix D: Local limit of f f -> S* amplitude 

In the local model with constituent quark masses m, ,the tri- 
angle quark-loop diagrams, depicted in Figs. 6b-c, reduce to 
the following expression 

s ioc-i(P 2 ^ 2 i^2) = m i8^ v h( m f)+ A ioc:i(p 2 ;q 2 uql) * 

xP^iquq^+B^^-qj^JP^iquq^, (D.18) 
where Ig(mf) is a gauge non-invariant term (constant) 



1 , 
. 2 . If, [ , mf-X 
Igyfitf) = - — - I dx\ I ax 2 



2n 2 



i 2 -X An 2 '' 



(D.19) 





X =X\{\ — X\ -x 2 )q\ +x 2 (l —x\ -x 2 )q\ + x\x 2 p 2 1 

which should be eliminated by suitable regularization, e.g., 
the Pauli-Villars regularization Ig(mf) — 7 g (Ap V ) = 0, and 



11 Bj (p 2 ;c/j is divergent in this kinematic. 



10 



the form factors read 



Aiocj (p 2 ;qi,ql) 



4(<71<?2)rc 



1 1— AT] 

, dx\ I dx 2 x 





4X-p 2 + q 2 {l-2x 1 )+q 2 (l-2x 2 ) 



mf-X 



(D.20) 



Bi oc j{p ;q\,q 2 ) 



l-x, 

dxi / dx 



(qm)n 2 J "" M J 



x 2 (l-2x 2 ) 



qmt-x) 



Note that, if one takes the local limit of the nonlocal ex- 
pression Eq. (C.15) by setting A^^ the contribution of 
nonlocal diagrams completely cancel a gauge noninvariant 
term. 

For special kinematics considered above, the form-factors 
become (1=1— x) 



Aioc-i(p ;p ,0) 



/' 2m? — p 2 xx ( 1 + 4xx) 
dx- 



\2% 2 



bLAp 2 >p 2 >°) = -^2 I dx 



(m? — xxp 
1 — 6xx 



2)2 



Jo mj — p 2 xx 



(D.21) 



Appendix E: Tensor structures for LbL amplitude 

Functions averaged over muon momenta can be represented 

as 



TsF=£(A),Z?7 



(E.22) 



where (A) ; - are the averages of scalar products with muon 
momentum in the numerator and muon propagators in the 
denominator (D! = (P + Qi) 2 + m 2 , D 2 = (P -Q 2 ) 2 +m 2 ) 



PQ 



2m 2 ! 



D 2 / 2m 2 



(A) 3 = { ) = (G1G2) 



-(Q1Q2) 



8m2 ' 

(l-^2) 2 

8,«2 



(E.23) 



(A>, 



1 



DjDz/ iHg|Gi||fi 2 |* 



arctan 



1-zf 



(A) 6 = <1) = 1, 

is the muon mass (P 2 = —m 2 ^ and 



t = 



(G1G2) 
IG1IIG2I 



4m 2 ! 

gT 



*=*|(1 -*)(!-*). 



(E.24) 



Zjy are polynomials in photon momenta 



V AA _ 
z l,l — 



J l,3 



(GiG 2 )(G 2 + (GiG 2 )), z^ = er 



3Gt + G? + G^ 



-Gi, z^ = ei 



»aa (QiQi)-Q\ 



-1.6 



Zff = (g(2m£ - Qf- (GiG 2 ))(Gf + (G1G2)), (E.25) 



7 AA _ 
z 2,l — 

V AA _ 
z 2,3 — 

7 AA 



Q 2 M-2Q 2 ) r7AA _Ql(Q 2 2 -2Ql) 



-1,2 



rt2 7 AA _ r,2 V AA _ r,2 Qi + Qi 
"kill z 2.4 —Uli ^2,6 — ~ J 



,AB 
J l,l 

7 AB 
"1,2 

,AB 
"1,3 
7 AB 



G?Gi((GiG 2 ) + G 



G3-G1 



Q 2 (6iG 2 ) 2 (G^ + (GiG 2 )) + GfG^ 



-G?(GiG 2 )((GiG 2 ) + Gl), 
z^4 =G2(GiG 2 )((GiG 2 ) + G2), 



,AB 
J l,5 



zt B 6 =((Q l Q2) z + QiQi 



2 ,(GiG2)±G| 
f (g^(G1 + (GiG 2 )(G5-G 2 ))- 



7 AB _ 
^2,1 — 



+ G?((GiG 2 )((GiG 2 ) + 5Gi) + 2G2: 



Z2 A f = #(G?(G^ + (GiG 2 )(G5-G^)H 



+ G2((GiG 2 )((GiG 2 ) + 5e?) + 2et) 

z 2 A f = G?((GiG2) + G?)((GiG 2 ) + G^ 
z 2 A ? = -G2((GiG2) + Gi)((GiG 2 ) + G2), 
7 ab _ n 2 n 2 n i (G1G2KG2 +G?) + 2G?G5 



AB _ (GlG2~(GlG2) XG2 + G 



2 

2 , ^2N 



-2.6 



G?G2Gi 



(E.26) 



(E.27) 
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